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INTRODUCTION 


Stirring a liquid doesn’t produce a flat surface but instead forms a dip in the middle and stirring it faster causes a 


greater depth. An adapted problem was proposed in the form of Newton’s Bucket as Newton mentioned:' 


“..the surface of the water will at first be flat, before it begins to recede from the centre and ascend up the 
sides of the bucket, forming itself into a concave figure, and the swifter the motion becomes, the higher will the 


DIL 


water rise, till at last, it becomes relatively at rest in it. 


The lack of a defined mathematical relationship between “swifter the motion” and “higher will the water 
rise” inspired an adapted experiment to quantify this observation and determine the relation between the 


synonymous angular velocity of the rotating fluid and the parabolic depth formed upon rotation. 





This topic is worthy of investigation because of its applications in our day-to-day lives. Rotating furnaces 
require a spinning liquid to create precise parabolic surfaces. These surfaces are utilised in fuel tanks of spinning 
rockets, solar cookers and centrifugal filters.““'” Astronomical telescopes make use of rotating fluids as focal length 
adjustable parabolic mirrors. The determination of an appropriate angular velocity is necessary to create parabolic 


surfaces with a precise depth.” 
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Figure 1 — image from the experiment which highlights Rigid Body Rotation. y-axis is 
the axis of rotation. w is the angular velocity of rotation. 





Theory 


Figure | is an example of rigid body rotation, where a substance rotates around a fixed imaginary axis of rotation — y-axis — 
consisting of all points inside the rotating body which remain stationary at a given angular velocity.“ The entire fluid 
moves as a rigid body such that individual particles are not deforming and have no shear stresses.“ The rotating fluid is 
stationary with respect to the rotating beaker and there is no relative motion between fluid particles; it appears that the 


water is at rest. 


It is assumed that the fluid at the beaker wall is moving with the same angular velocity as the beaker and adhesion 
of liquid to the beaker is assumed large enough to neglect any slipping at the wall. However, in reality, there are frictional 
forces present between the liquid and the surface of the rotating vessel. There exists a boundary layer near the rotating 


vessel where the fluid velocity in that region is different from the rest of the fluid. “# 


The derivation of the relationship between angular velocity of the fluid and its parabolic depth is considered in an 
inertial frame of reference since the consideration of a rotating (non-inertial) frame of reference along with hydrostatic 


29 1x 


force can produce “paradoxical results”. 


The following derivation assumes the consideration of Newtonian fluids only since Non-Newtonian fluids require 
a yield stress to overcome before rotational motion is initiated, which is not the case for Newtonian fluids. Additionally, 


resistance to free flowing motion (viscosity) is assumed to be sufficiently low.** 


The high number of assumptions in the model is a limitation of the experiment as not all of them may be met 


while conducting the experiment. 
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Derivation of the Relationship 


Figure 2— Force Body Diagram 





Figure 3 — Magnified Force Body Diagram 


To determine forces on the rotating fluid, the forces a single droplet of water (denoted by the black circle) is 


considered: 
The tangent in Figure 2 and 3 is used a reference line to mark angles. 


Primarily, there are just two forces acting on the droplet of water in circular motion; first is the weight of the 


liquid acting downwards: 


Fy = mg () 
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m is the mass of the water droplet and g is the acceleration due to gravity. 


The second force is the force created due to fluid pressure — N — which acts perpendicular to the tangent as seen in 
Figure 2 and 3. Inside the fluid, the pressure force is exerted isotropically — equally in all directions. However, this force 
becomes unbalanced at the parabolic surface since the force can only push the particles from inside the fluid. Hence, at the 


parabolic surface, the force created due to fluid pressure is directed outward and perpendicular to the surface.*# 


A centripetal force (F.) must act on the water droplet in circular motion with its centre on the y-axis. F, is not 


another force but is a component of N acting towards the centre: 


F. = Nsin(6) (2) 


Nsin(®) is unbalanced (no opposing force) so it alone provides centripetal acceleration to keep the liquid moving 


in a circular path with angular velocity . The liquid accelerates and F, is also given by: 


F, = mw*Ar (3) 


Ar is the distance from the centre of the circular path to the droplet of water. 
Equating (2) and (3): 


N sin(@) = mw?r (4) 


After the shape is formed, the weight of the particle, (1), and the cosine component of N, Ncos®, are equal in 


magnitude but opposite in direction because the water droplet doesn’t undergo any vertical displacement: 








N cos(@) = mg (5) 
Dividing (4) with (5): 
N sin@ mw?r 
NcosO mg (6) 
wr 
tan(@) = ae (7) 
The right angled triangle in Figure 3 suggests: 
Ah _ dh 
tan(@) = oo ae (8) 


Equating (7) and (8) gives the differential equation: 
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dh wr 
dr. g ) 
dic ay 

g 


Integrating both sides: 


Jah =f" ar 


h and r are integrated: 


fadh==frar 


The general solution to this differential equation is: 





(10) 


h is the vertical displacement of an arbitrary point from the base of the container. 


ho is the constant of integration and is the height from the base of the container to the minimum of the parabola 
formed. 


In context of the research question, (10) is rearranged: 





7 (1) 


The parabolic depth, hg is the vertical height difference between the minimum point of the parabola and the 
height of the surface at r=ro and is: *“ 


hg=h—-hy (12) 
www.Gpre.org 
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Parabolic Depth = 


hg =h-hy 


Mm ———— 
r=" 


Figure 4 — highlighting all the heights that will be 
discussed 





The different heights to be measured are: 
h is the vertical displacement for the point at r=ro. 
The conditions of this experiment were: 
e hg < hand 
e hj)=O0 
These conditions imply that the minimum of the parabola must be above or at the surface of the beaker. 


(12) is substituted in (11) andr is substituted with ro— radius of the beaker: 





ha = (13) 
Relating the angular velocity (@) to the parabolic depth (ha): 
hg = 2-0? 14 
tog (14) 
Replacing the constants (ro and g): 
hg = kw? (15) 
hq « w? (16) 
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The relationship to be investigated is shown in the equation of a parabola which relates angular velocity (w) of the 
fluid to the parabolic depth (hg). It is predicted that plotting ha against ? should ideally give a straight line passing through 


the origin with a gradient given by the coefficient of w” from (14). 


The reasoning behind the direct relationship between w* and hg, from the Research Question, is understood 


through application of pressure: 


p= ae (17) 


p is pressure on fluid particles and p is the density of the fluid.*“"™ 
This equation also demonstrates the parabolic nature of the curve. It shows that the pressure for points at the same 


distance, r, away from the axis of rotation, is the same. Thereby, the fluid rotating with higher angular velocity will have a 


higher pressure than the fluid rotating at a lower angular velocity (Ap is positive). 


Now, the equation above is re-arranged in terms of the variables defined to relate the fluid pressure to the 


parabolic depth: 





p= pg x 


wr? 
29 


p= pghq (18) 
Since g is a constant and the density of the liquid remains same: 
pxhg 


The pressure exerted on fluid particles is directly proportional to the parabolic depth. Thereby, (17) showed that 
the increasing angular velocity increases the pressure exerted on the fluid particles which causes a greater parabolic depth 


(ha), seen in (18). 
Experiment 


The aim of the experiment was to change the angular velocity (@) of the fluid (which was done by changing the input 


voltage to the motor) and determine the parabolic depth (ha) measured by fixing a ruler next to rotating beaker. 


Controlled variables are present to ensure that the angular velocity is the only factor which changes parabolic 
depth. These included: radius (r) of the beaker, volume of water used, the time after which the height levels were 


measured, and the general setup. 
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Retort stand <—— 
White tile 


Clamp 


Beaker Ruler 


containing 
Wooden stand mounted 
on the motor’s axle 


Power 
Supply 


Black Crocodile 
Clip wire Red Crocodile 
Clip wire 
Laptop for 
Stopwatch data collection 


Phone for 
video recording 


Figure 5 — Diagrammatical image of the setup 





The diagram of the setup is seen in figure 5. 


The basic choice of materials from the school laboratory ensures the experiment can be repeated with ease. 


Figure 6 — image of the actual setup 
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The actual setup is seen in figure 6. 


—_ ——_—— White Tile 


Beaker 


Adhesive 
Wooden disk 


2 washers 


Figure 7 — Close view of the motor, washers, wooden 
disk, adhesive, and beaker 





Apparatus List 


5 retort stands with 5 clamps — to fix apparatus as seen in figure 6. Circular wooden disk (radius 5.3 cm + 0.1 cm) as the 
base of the beaker attached to it using an adhesive. 2 washers were used to fill the gap of the axle. Beaker containing 100 
ml of water. Power supply used to alter the voltage supplied to the motor which changes the angular velocity. Stopwatch 
used while calculating angular velocity of the rotating beaker. (+0.01s). Motor — to rotate the wooden disk and the 


connected beaker. A close up view of the rotating beaker is shown in figure 7. 
Procedure 


The range of independent variables had to be decided which followed a criterion. The lowest angular velocity (8.2 rad s“!) 
should produce a parabola with non-negligible depth such that measurement is possible. The highest angular velocity (33.3 
rad s') should produce a complete parabola such that its minimum is above or at the base of the beaker. The other angular 
velocities must be between these limits and its overall range (in rad s‘!) is: 8.2, 11.9, 17.6, 22.4, 30.2, 33.3, and 39.6. 
However, one selected angular velocity (of 39.6 rad s! whose data was discounted) was above the upper limit to evaluate 


its parabolic shape. 


Before determining the angular velocity, the time period was calculated using: 


. F 10 seconds 
Time Period = T = ————_ (19) 
Revolutions in 10 seconds 
@ was calculated: 
21 


Here, it was earlier assumed that the angular velocity of the rotating beaker is the same as that of the rotating 
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water. However, deviations from this assumption limit the accurate determination of the angular velocity values. 


The voltage was first set such that the angular velocity was 8.2 rad s!. Time of 1 minute was given for the water 
to come at rest with respect to the rotating beaker. The heights (ho and h from figure 4) were measured and the parabolic 


depth (ha) was calculated. 


Lee in 10 lL ! hie (col) 0.05 om h (cm) + 0.05 cm 
seconds 
—_s s)+ Height at rest 
as Sal paooo sins eee 





























Table 1 - Raw Data 


Data Collection 





hg (cm) + 0.1 cm 





angular 
Time Period | Velocity @ 
(s)+0.01s | (rads) + 

0.01 rad s" 





Sr No. 





09109108 
1.6] 1.6| 1.6 


0.19 % a6 1100 re , 2 re a a rT 
0.16 1600 99 NA|N/AIN/AIN/AIN/A N/A NA 


Table 2 - es Data calculated using data from Table 1 




















SIDA a ee lb le 





Table 1 show cases all data that was measured and further calculations are based on this data only. 


Reasons for “N/A” and yellow highlighted cells are explored later. Time period, angular velocity and ha were 


calculated using (19), (20) and (12) respectively. 
Data Analysis 


Firstly, o = 8.2 rad s! had the greatest value of error in hg and is greater than the instrumental error. This was because the 
parabola that was formed had negligible depth, making its measurement with a ruler difficult as the ruler’s least count is 


0.1cm and doesn’t allow measurement of heights to great accuracy. 


@ = 33 rad s! had its minimum of parabola located at the base of the beaker, quantitatively seen by the hoaverage 


being 0.0 cm (rounded to 1 decimal place) and also seen in figure 8. 
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Figure 8 — image of the parabolic surface at 33 rad s-1 when the minimum 
of the parabola is on the surface of the beaker, making this angular 
velocity the upper limit of the range of angular velocities selected 


ee 


of 


‘hy>h and 
h)<0.0 cm 


Figure 9 — image of the parabolic surface for an 
angular velocity of 40 rad s-1 whose hypothetical 
minimum would lie below the surface of the beaker 
if the parabola is extended. 





As mentioned already, @ = 40 rad s! was chosen with a value greater than the maximum angular velocity, with a 


shape shown in Figure 9. 


The parabola was extended to visualise a minimum point. Even though the depth of the extended parabola could 
be determined using image editing software, it would be a non-viable option as it surpasses the conditions of the 
experiment (hg < h and hy = 0). This explains ‘N/A’ in Table 1 for ho and ha. Thereby, the range of angular velocities a 


limitation of the experiment which can be improved by using a beaker of the same radius but of a larger height such that 
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parabolas of greater depths can fit into it, allowing a higher range of independent variables. 


The two yellow highlighted cells in Table 1 are two anomalies of the experiment because their values deviate 
compared to other values of ho and h. Thereby, these values were discounted while calculating hg, explaining ‘N/A’ in hal 
and hg2 for the w? = 910 rad? s? in Table 3. This anomaly could exist if the heights were measured before the rotating fluid 


became stationary with respect to the beaker. 


Error in ha for @ (in rad s“!) values of 18 and 22 were low due to the non-negligible parabolic depth that could be 


measured easily because of low turbulent motion. 


As «0? (in rad’ s*) increases from 67 towards 1100, hg (in cm) also increases from 0.6 towards 6.2 (from Table 3). 


The positive correlation between w” and hy is as expected from the prediction. 


Straight Line Graph Analysis 


Graph 1 - Parabolic Depth (hy) versus Angular Velocity 


Squared (w?) 
7.00 
h, = 0.0057w? + 0.0929 


6.00 = SS 


5.00 ve 
hy = 0.0051? + 0.1350 








ot atan ase cef® R? = 0.9842 
= 4.00 Ros 
a TT 
= 3.00 hy =0.0051w? + 0.4317 
2.00 
1.00 5 a 
0.00 
0 200 400 160% .-2, 800 000 200 
. ue am w2 (RP 52) 808 anes ete 


Plotting w? versus hg: 


The vertical error bar was of a total height of 0.30 cm. The highest error in hg (0.15 cm for w@=8.2 rad s! ) was 


applied to all the other data points as the other error values were too small to be shown. 


As predicted, the graph of w? versus hg produces a straight line. However, the line doesn’t pass through the origin 
but instead passes through a y-intercept value of 0.135. Even though systematic errors exist (=0.135) they are negligible 
and a better methodology will improve the accuracy of the results. A value of R’ close to 1 shows a strong correlation 


between the two variables and implies that most of the data points can be approximated by the trend line that was created. 
It was predicted that the value of the gradient should be: 


2 

a rT 
gradient = 
29 
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The theoretical value is calculated by replacing ro=3.25 cm and g=980.67 cm s” in the above equation: 


3.252 
2980.67 





gradient = Mtheoretical = 


Meheoretical = 9.005385 ~ 0.0054 


The experimental value of the gradient is: 
max gradient = 0.0057 
min gradient = 0.0051 


0.0057 — 0.0051 


error in gradient = ——_— 0.0003 


Mexperimental = 0.0053 + 0.0003 


The closeness of the experimental and theoretical values suggests the high accuracy of the data received, implying 


a strong direct relationship between w* and ha and showing that no other variables significantly affected ha. 


Logarithmic Analysis 


A logarithmic analysis verifies the power to which the angular velocity is related to the parabolic depth. The general 


positive correlation assumes: 
n 
hg X@ 


ha = Kw" 


Taking the natural log on both sides: 


Inhg = In(Kw") 


Using the properties of logarithms: 


In(hg) = nin(@w) + In(k) (21) 


This is synonymous to the equation of a straight line: 
y=Mxt+c 


Thereby, plotting In(ha) against In(@) should ideally have n=2 and the value of In(K) shows the contribution of the 


controlled variables. 
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@ (rad s') +|ha average (cm) : } error in 

Sr No. itimde LO teh In(@)| In(ha) |hamax| hamin In(ha) 

1 8.2 0.6 2.10|-0.48 | 0.8 | 0.5 0.24 

2 12 0.9 2.48|-0.11}| 1.0 | 0.8 0.11 

3 18 1.6 2.87| 0.44 | 1.6 | 1.5 0.03 

4 23 2.8 3.12| 1.03 | 2.8 | 2.8 0.00 

5 30 44 3.41| 148] 45 | 43 0.02 

6 33 6.2 3.51 | 1.82 | 6.2 | 6.1 0.01 

7 40 N/A 3.68| NA | N/A | N/A NA 
































Table 4: Data for logarithmic analysis 


Here again, “N/A” is written for @ = 40 rad s”' as its ha values are not considered. 


In hgmax-In hgmin 
2 


error in In(hg) = 
Sample Calculation 


Sample calculation for m = 8.2 rads": 


In(0.8)—In(0.5) _ 


error in In(hg) = 0.23500 = 0.24 


Similar calculations were conducted for other values of o. 


Graph 2 - In(hg) versus In(w) 





In(hy) = 1.968In(w) - 4.847 


2.50 
Infhy) =1.635In(w) - 4.065 

2.00 “ft R2 = 0.9783 
1.50 aye eal 

all In{hy)=1.299In(w) -2.972 
+00 - 
0-50 et 
0.00 _ ai aa 

“so Poa 

2.00 at 50 3.00 3.50 
-0.50 ya 
-1.00 


Plotting In(ha) against In(w): 


The equation of the best fit line is: 


In(w) 


In(hg) = 1.635 In(w) — 4.065 


The gradient is 1.64 when rounded to 3 significant figures. 
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max gradient = 1.968 


min gradient = 1.299 


1.968-1.299 


error in gradient = ———s 0.3345 = 0.335 


average gradient = 1.635 + 0.335 


The upper bound for the value accurately predicts the quadratic relationship. The gradient is 1.64 + 0.33 when 


rounded to 2 decimal places and this estimates the value of n in (23). The y-intercept is: 
c = In(K) = —4.065 


Kee? = 0.117 


0.017 is the constant which represents the factor that the controlled variables contribute to influence the parabolic 


depth. 
The logarithmic analysis suggests: hg = 0.017w*** which is an 18% deviation from w?. 
Evaluation 


Some of the reasons which can explain the deviation above are listed in this section, suggested improvements are also 


listed. 


The human reaction error influences the value of time period to determine the angular velocity. This was tackled 
by video recording the stopwatch along with the rotating beaker. This allowed the accurate determination of the number of 


rotations in a 10.00s period. This can further be improvised by using a photocell connected to a computer. 


The point in time at which the height measurements were taken is subjective as enough time had to be given for 
the rotating water to be stationary with respect to the rotating beaker. The subjectivity in the determination of this point 
includes inevitable errors. And the early measurement of the depth can lead to misleading values (cells highlighted in 


yellow in Table 1). 


Additionally, difficulties to measure the heights increased when motion got turbulent at higher angular velocities 
as the beaker began to wobble. This causes a random error seen at w* = 910 rad? s*, which is away from the trend line in 


Graph 1, this reduces the precision of the results. 


The angular velocity was kept the same by keeping a constant input voltage. However, during the process of data 
collection, the input voltage, instantaneously, had minor fluctuations (for example from 1.1V to 1.2V and back to 1.1V). 
These fluctuations along with resistance in wires can cause negligible changes in angular velocity which directly impacts 


the parabolic depth. The suggested improvement can be to use a variable speed motor. 


A general improvement to correct the errors above is to take more repeats to reduce the errors’ effect on the precision of 


the results. 


The y-intercept in Graph | is the value of the systematic errors since these errors cause a fixed deviation of the 


results and cause the line to pass through an intercept rather than the origin (as expected in the prediction assuming a 
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perfect experiment). The general method to improve these errors is to utilise a more sophisticated setup that tends to be 


difficult, expensive and cumbersome to setup in a school laboratory. *¥##** 
CONCLUSIONS 


In conclusion, by means of this experiment, the relationship between the angular velocity of a fluid in rigid body rotation 
and its parabolic depth was experimentally determined. The theory suggested a linear relation between hg and w’, however, 
experimentally a value of w!°* was determined — a 18% deviation. This deviation is explained by the errors, limitations and 
assumptions in this experiment which inhibits maximum accuracy and precision of the results. The final relationship 
received through the logarithmic analysis was: hg=0.017@!'. This was explained through the direct relationship between 


angular velocity and pressure exerted on fluid and consequently a higher parabolic depth. 


The biggest factor which leads to this deviation is the basic choice of apparatus used which, as mentioned already, 
can be used in any school laboratory. Even though these two values (1.64 and 2) are similar, the difference can be 


minimized by avoiding a few errors to improve the accuracy and precision of the results achieved. 
Even though the RQ is answered, further experiments (amongst others) to study Rotating Fluids can be to study: 
e =©Them in a rotating (non-inertial frame) of reference. 


e Its applications in parabolic mirrors to determine its focal length. 
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